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ABSTRACT 

Using recent measurements of the supermassive black hole (SMBH) mass function, we find that SMBHs 
are the largest contributor to the entropy of the observable universe, contributing at least an order of magnitude 
more entropy than previously estimated. The total entropy of the observable universe is correspondingly higher, 
and is Sobs = 3-l!l;i 7 x 10'""^^. We calculate the entropy of the current cosmic event horizon to be Sqeh = 
2.6 ±0.3 X 10^'^\ dwarfing the entropy of its interior, Sceh int = l ^^a' x 10'°'' A;. We make the first tentative 
estimate of the entropy of weakly interacting massive particle dark matter within the observable universe, 
5'dm = 10^^^^ k. We highlight several caveats pertaining to these estimates and make recommendations for 
future work. 

Subject headings: black hole physics — cosmology: miscellaneous — diffusion — elementary particles — 
gravitation — neutrinos 



1. INTRODUCTION 

The entropy budget of the universe is important because 
its increase is associated with all irreversible processes, on 
all scales, across all facets of nature: gravitational clustering, 
accretion disks, supernovae, stellar fusion, terrestria l weather, 
and chemical, geological and bi ological processes (iFrautschil 
[TQSl lLineweaver & Eganll2008h. 

Rec ently, iFrampton et al.l (120091) and iFrampton & Kepharll 
(120081) reported the entropy budget of the observable uni- 
verse. Their budgets (listed aside others in Table [U esti- 
mate the total entropy of the observable universe to be Sohs ^ 
lO'^'A:- 10'""'^:, dominated by the entropy of supermassive 
black holes (SMBHs) at the centers of galaxies. That the in- 
crease of entropy has not yet been capped by some limiting 
value, such as the holographic bound (n; Hooft 1993; Sussking 
11995b at 5n,ax ~ 10'~^/fc ( IFrampton e t al. 2009), is the reason 
dissipative processes are ongoing and that life can exist. 

In this paper, we improve the entropy budget by using re- 
cent observational data and quantifying uncertainties. The pa- 
per is organized as follows. In what remains of the Introduc- 
tion, we describe two different schemes for quantifying the in- 
creasing entropy of the universe, and we comment on caveats 
involving the identification of gravitational entropy. Our main 
work is presented in Sections|2]and[3] where we calculate new 
entropy budgets within each of the two accounting schemes. 
We finish in Section|4]with a discussion touching on the time 
evolution of the budgets we have calculated, and ideas for fu- 
ture work. 

Throughout this paper we assume flatnes s {ilk = 0) as pre- 
dicted by infl ation (iGuthlll 98 lULindd 19821) and supported by 
observations (ISper gel et al .1 120071) . Adopted values for other 
cosmological parameters are /z = 0.705 ± 0.013, ut, = ilth~ = 
0.0224 ± 0.0007, uj,„ = flji^ = 0.136 ±0.003 (iSehak et al.1 



l2006h . and Tcmb = 2.725 ± 0.002 K (iMatheret al.1 [19991: 
quoted uncertainties are la). 

1.1. Two Schemes for Quantifying the Increasing Entropy of 
the Universe 

Modulo statistical fluctuations, the generalized second law 
of thermodynamics holds that the entropy of the universe (in- 
cluding Bekenstein-Hawking entropy in the case of any re- 
gion hidden behind an event horizon), must not decrease with 
time ([Bekenstein 1974; Gibbons & Hawking 1977). Within 
the FRW framework, the generalized second law can be ap- 
plied in at least two obvious ways: 

1 . The total entropy in a sufficiently large comoving vol- 
ume of the universe does not decrease with cosmic 
time, 

^'^comoving volume 0- (1) 

2. The total entropy of matter contained within the cosmic 
event horizon (CEH) plus the entropy of the CEH itself, 
does not decrease with cosmic time. 



interior 



(2) 
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In the first of these schemes, the system is bounded by a 
closed comoving surface. The system is effectively isolated 
because large-scale homogeneity and isotropy imply no net 
flows of entropy into or out of the comoving volume. The 
time-slicing in this scheme is along surfaces of constant cos- 
mic time. Event horizons of black holes are used to quan- 
tify the entropy of black holes, however the CEH is neglected 
since the assumption of large-scale homogeneity makes it 
possible for us to keep track of the entropy of matter beyond it. 
A reasonable choice for the comoving volume in this scheme 
is the comoving sphere that presently corresponds to the ob- 
servable universe, i.e., the gray area in Figure[T] Correspond- 
ingly, in Section|2]we calculate the present entropy budget of 
the observable universe and we do not include the CEH. 
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The second scheme is similar to the first in that we time- 
slice along surfaces of constant cosmic time. However, here 
the system (yellow shade in Figure [T]) is bounded by the 
time-dependent CEH instead of a comoving boundary. Mi- 
gration of matter ac ross the CEH is not negligible, and the 
CEH entropy (G ibbons & Hawking|ll977|) must be inc luded 
in the budget to account for this (e.g. lDavis etaLll2003h . The 
present entropy of the CEH and its interior is calculated in 
Section[3] 
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Fig. 1 . — These two panels show the pailicle horizon (see Equation |42]and 
Figure |9) and the cosmic event horizon (see Equation I46t as a function of 
time. The difference between the two panels is the spatial coordinate system 
used: the ,v-axis in the bottom panel is proper distance D and in the top panel 
it is comoving distance x = ^ ■ where a is the cosmic scalefactor. The ori- 
gin is chosen so that our galaxy is the central vertical dotted line. The other 
dotted Unes represent distant galaxies, which are approximately comoving 
and recede as the universe expands. The region inside the particle horizon 
is the observable universe. The comoving volume that corresponds to the 
observable universe today, about 13.7 Gyr after the big bang, is filled gray. 
In scheme 1 , the entropy within this comoving volume increases (or remains 
constant) with time. Alternatively, in scheme 2 the entropy within the event 
horizon (the region filled yellow), plus the entropy of the horizon itself, in- 
creases (or remains constant) with time. 



1.2. Entropy and Gravity 

It is widely appreciated that non-gravitating systems of par- 
ticles evolve toward homogenous temperature and density 
distributions. The corresponding increase in the volume of 
momentum-space and position-space occupied by the con- 
stituent particles represents an increase in entropy. On the 
other hand, strongly gravitating systems become increasingly 
lumpy. With "lumpyness" naively akin to "orderliness", it is 
not as easy to see that the total entropy increases. In these 
systems the entropy is shared among numerous components, 
all of which must be considered. 

For example, approximately coUisionless long-range grav- 
itational interactions between stars result in dynamical re- 
laxation of galaxies (whereby bulk motions are dissipated 
and entropy is transferred to stars in the outer regions of 



the galaxy; iLvnden-Belll 1 1 967l) and stellar evaporation from 
galaxies (whereby stars are ejected altogether, carrying with 
them energy, angular momentum and ent ropy, and allow- 
ing w hat remains behind to contract; e.g. Bi nney & Tremaina 
l2008h . In more highly dissipative systems, i.e., accretion 
disks, non-gravitational interactions (viscosity an d/or magne- 
torotational instability; iBalbus & HawlevI 120021) transfer an- 
gular momentum and dissipate energy and entropy. 

In addition to these considerations, entropy also increases 
when gravitons are produced. A good example is the in-spiral 
of cl ose binaries, such as the Hulse-Taylor binary pu lsar sys- 
tem (iHulse & TayloiilI975l:rWeisberg & Tayloiil2005h . Gravi- 
tational waves emitted from the system extract orbital energy 
(and therefore entropy) allowing the system to contract. 

The e ntropy o f a ge neral gravita tional field is still not 
known. iPenrosd (1 1 9871 119791 12004 has proposed that it is 
related to the Weyl curvature tensor W^^nX- In conformally 
flat spacetimes (such as an ideal FRW universe), the Weyl 
curvature vanishes and gravitational entropy is postulated to 
vanish (to limits imposed by quantum uncertainty). In clumpy 
spacetimes the Weyl curvature takes large values and the grav- 
itational entropy is high. While Ricci curvature 7?^^ vanishes 
in the absence of matter, Weyl curvature may still be non-zero 
(e.g. gravitational waves traveling though empty space) and 
the corresponding gravitational entropy may be non-zero. 

If these ideas are correct then the low gravitational en- 
tropy of the early universe comes from small primordial grav- 
itational perturbations. Gravitational entropy then increases 
with the growing amplitude of linear density fluctuations pa- 
rameterized through the matter power spectrum P{k). The 
present gravitational entropy, however, is expected to be dom- 
inated by the nonlinear overdensities (with large Weyl ten- 
sors) which have formed since matter-radiation equality. 

In extreme cases, gravitational clumping leads to the forma- 
tion of black holes. The ent ropy of black holes is well know n 
(iBekensteinll 1971 iHawkin g 1976: Strominger & Vafal[l996l) . 
The entropy of a Schwarzschild black hole is given by 



kc^ A 47rkG 



on 4 



(3) 



where A -- 



\6tiC-M- 



is the event-horizon area and M is the black 



hole mass. 

Because gravitational entropy is difficult to quantify, we 
only include it in the two extremes: the thermal distribution 
of gravitons and black holes. 

2. THE PRESENT ENTROPY OF THE OBSERVABLE UNIVERSE 

The present entropy budget of the observab l e uni verse 
was estima ted most r ecently by iFrampton et al.l ( 120091) and 
lF?ampton & KephartI (l2008h. Those pap e rs and earlier 
work dKolb & Turned n^SlTlFrautschil [T98a IPenrosd 120041: 
iBousso et al.ll2007f) identified the largest contributors to the 
entropy of the observable universe as black holes, followed 
distantly by the cosmic microwave background (CMB) and 
the neutrino background. The last column of Table[T]contains 
previous estimates of the entropy in black holes, the CMB and 
neutrinos, as w e ll as several less significant components. 

Sections r2.1| - |2.7| below describe the data and assumptions 
used to calculate our entropy densities (given in Column 2 
of Table [Til. Our entropy budget for the observable universe 
(Column 3 of Table [Til is then found by multiplying the en- 
tropy density by the volume of the observable universe Vobs, 



(4) 
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where i, is the entropy density of component /. The volume 
of the observable universe is (see Appendix) 



Vobs = 43.2 ± 1.2 X 10^ Glyr^ 
= 3.65 ±0.10 X lO'^^OT^ 



(5) 



2.1. Baryons 

For a non-relativistic, non-degenerate gas the specific en- 
tropy (entrop y per baryon ) is given by the Sakur- Tetrode 
equation (e.g. lBasu & Lynden-Bell 1990i) 



is /lib) = — y]n,ln 



ZiiT)i2TTmikT)ie--nJ^h~ 



(6) 



where ; indexes particle types in the gas, n, is the j'th particle 
type's number density, and Z,(r) is its internal partition func- 
tion. IBasu & Lynden-Bell (1990) found specific entropies be- 
tween 11 ^: and 21 k per baryon for main-sequence stars of 
approximately solar mass. For components of the interstellar 
medium (ISM) and intergalactic medium (IGM) they found 
specific entropies between 20 k {Ho in the ISM) and 143 k 
(ionized hydrogen in the IGM) per baryon. 

The cosmic entropy density in stars can be estimated by 
multiplying the specific entropy of stellar material by the cos- 
mic number density of baryons in stars «/„ : 



m„ 



3H^ rj. 



SttG m„ 



(7) 



Using the stel lar cosmic de nsity parameter fl^ = 0.0027 ± 
0.0005 (Fuku gita & PeeblesI 12004). and the range of spe- 
cific entropies for main-sequence stars around the solar mass 
(which dominate stellar mass), we find 



5,=0.26±0.12yt;7r', 
5* =9.5 ±4.5 X 10^° k. 



(8) 
(9) 



Similarly, the combined energy density for the ISM and IGM 
is f7gas = 0.040 ±0.003 (iFukugita & PeeblesI l200l . and by 
using the range of specific entropies for ISM and IGM com- 
ponents, we find 



Sgas = 20± 15 k m ^, 
5gas=7.1±5.6x lO^'/t. 



(10) 
(11) 

The uncertainties in Equations (|9]l and ( fTTT i are dominated by 
uncertainties in the mass weighting of the specific entropies, 
but also include uncertainties in il*, ilgas and the volume of 
the observable universe. 

2.2. Photons 

The CMB photons are the most significant non-black 
hole contributors to the entropy of the observ able universe. 
The distribution of CMB photons is thermal dMather et alj 
119941) with a prese nt temperature of = 2.725 ± 0.002 K 
(iMather etal^l999^ . 

The entropy of the CMB is calculated u sing the equation 
for a black body (e.g. iKolb & Turned (I1990I) ). 



It:' k" 



45 c^r?''^"' 
= 1.478 ±0.003 X 10'*/tm-^ 
5^ = 2.03±0.15 X lO^H, 



(12) 



(13) 



where = 2 is the number of photon spin states. The uncer- 
tainty in Equation (fT3T l is dominated by uncertainty in the size 
of the observable universe. 

The non-CMB photon contribution to the entropy budget 
(including starlight and heat emitted by the ISM) is somewhat 
less, at around lO^^ k jFrautschil Il982t iBousso et all 120071; 
iFrampton et"ani2009h . 

2.3. Relic Neutrinos 

The neutrino entropy cannot be calculated directly since the 
temperature of cosmic neutrinos has not been measured . Stan- 
dard treaties of the radiation era (e.g. IKolb &Turnei1 [19901: 
lPeacockll999h describe how the present temperature (and en- 
tropy) of massless relic neutrinos can be calculated from the 
well known CMB photon temperature. Since this background 
physics is required for Sections [2.41 and |231 we summarize it 
briefly here. 

A simplifying feature of the radiation era (at least at known 
energies < lO'^eV) is that the radiation fluid evolves adiabati- 
cally: the entropy density decreases as the cube of the increas- 
ing scalefactor i^d oc a~^. The evolution is adiabatic because 
reaction rates in the fluid are faster than the expansion rate H 
of the universe. It is convenient to write the entropy density 
as 



27r2 k^ 3 _3 



(14) 



45 cHi^' 

where g^s is the number of relativistic degrees of freedom in 
the fluid (with m < kT / c~) given approximately by 

bosons, i 



fermions. J 



For photons alone, g^,s = = 2, and thus Equation (O be- 
comes Equation (T% . For photons coupled to an electron- 
positron component, such as existed before electron-positron 

annihilation, g*s = g-, + lge± = 2+l4= y- 

As the universe expands, massive particles annihilate, heat- 
ing the remaining fluid. The effect on the photon temperature 
is quantified by inverting Equation fT4[ . 



-1/3 



(16) 



The photon temperature decreases less quickly than a"' be- 
cause g^s decreases with time. Before electron-positron e^ 
annihilation the temperature of the photons was the same as 
that of the almost completely decoupled neutrinos. After e^ 
annihilation, heats only the photons, the two temperatures dif- 
fer by a factor C, 



(17) 



A reasonable approximation C « (4/1 1)'/-' is derived by as- 
suming that only photons were heated during annihilation, 
where 4/11 is the ratio of g^,s for photons to g^,s for photons, 
electrons, and positrons. 

Corrections are necessary at the 10"^ level because neu- 
trinos had not completely decoupled at annihilation 
jGnedin & Gnedi n 1998). The neutrino entropy density 
is computed assuming a thermal distribution with = 
(4/1 l)'/-'r.y, and we assign a 1% uncertainty. 



2_TrJ^ 
45 cHi^ 
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TABLE 1 

Current Entropy of the Observable Universe (Scheme 1 Entropy Budget) 



Component 



Entropy Density i [k m 



Entropy S [k] Entropy S [k] (Previous Work) 



SMBHs 


8.4+15 ^ 10^^ 


3.1+3 « X lO""* 


10""[1], lO'^'P], lO'^^P] 


Stellar BHs (2.5-15 Mq) 


1.6xl0"^:2 


5.9x10"-?:! 


lo'^'m, 10'** [4] 


Photons 


1.478 ±0.003 X 10' 


5.40±0.15 X 10**' 


10'*'*[L2,4], 10'*'[5] 


Relic Neutrinos 


1.411 ±0.014 X 10' 


5.16±0.15 X 10'*' 


10*"*[2],10'*'[5] 


WIMP Dark Matter 


5 X 10'±' 


2x 1088±' 




Relic Gravitons 


1.7 X 10^-2.5 


6.2 x 108^-?:5 


10'*'' [2, 3] 


ISM and IGM 


20 ±15 


7.1 ±5.6 X 10'*' 




Stars 


0.26 ±0.12 


9.5 ±4.5 X 1(P 


10™ [2] 


Total 


S.4+l^ X 10''^ 


X 10™ 


10""[1], lO^'-'P], lO'^^P] 


Tentative Components: 
Massive Halo BHs (10^ Mq) 


1025 


10106 


lO'^M 


Stellar BHs (42-140Mq) 


,0+0.8 

8.5 X 10'**-i 6 


3.1 X 10"-?:' 





Note. — Our budget is consistent with previous estimates from the literature with the exception that SMBHs, 
which dominate the budget, contain at least an order of magnitude more entropy as previously estimated, due to 
the contributions of black holes 100 times larger than those considered in previous budgets. Uncertainty in the 
volume of the observable universe (see Appendix) has been included in the quoted uncertainties. Massive halo 
black holes at lO' Mq and stellar black holes in the range 42—140 Mq are incl uded tentatively sinc e their existence 
is spec ulat ive. They are not counted in the budget totals. P revi ous work: [11IPenros3 )2004l), [21 IFrampton et alj 
f2009l) . 131 IFrampton & Kepharll {2003) . 141 IFrautschil (190) . 151 IKolb & TumeiH1981h . 161 IFramptoij I2009bl) . 



= 1.411 ±0.014 X lOHm- 



(18) 



where gu = 6{3 flavors, 2 spin states each). The total neutrino 
entropy in the observable universe is then 



5„=5.16±0.14x 10'*'' ;t 



(19) 



with an uncertainty dominated by uncertainty in the volume 
of the observable universe. 

Neutrino oscillation experiments have demonstrated that 
neutrinos are massive by measurin g differences between 
the three n eutrino mass eigenstates (C leveland et alJ 119981 : 
lAdamson et a l. 2008; Abe et al. 2008). At least two of the 
mass eigenstates are heavier than ^ 0.009 eV. Since this 
is heavier than their current relativistic energy (^ C Tj = 
0.0001 eV; computed under the assumption that they are 
massless) at least two of the three masses are presently non- 
relativistic. 

Expansion causes non-relativistic species to cool as in- 
stead of a~\ which would result in a lower temperature for the 
neutrino background than suggested by Equation (fTTT i. The 
entropy density (calculated in Equation [TSll and entropy (cal- 
culated in Equation \T% are unaffected by the transition to 
non-relativistic cooling since the cosmic expansion of rela- 
tivistic and non-relativistic gases are both adiabatic processes 
(the comoving entropy is conserved, so in either case s oc a~^). 

We neglect a possible increase in neutrino entropy due to 
their infall into gravitational potentials during structure for- 
mation. If large, this will need to be considered in future 
work. 

2.4. Relic Gravitons 

A thermal background of gravitons is expected to exist, 
which decoupled from the photon bath around the Planck 
time, and has been cooling as T^^^ oc since then. The 
photons cooled less quickly because they were heated by the 
annihilation of heavy particle species (EquationfTSTl. Thus we 
can relate the current graviton temperature to the current pho- 
ton temperature 



T = 

i gray — 



^*s(fplanck) 



1/3 



(20) 



where g*s(fpianck) is the number of relativistic degrees of free- 
dom at the Planck time and g^sito) = 3.91 today (this is ap- 
propriate even in the case of massive neutrinos because they 
decoupled from the photon bath while they were still relativis- 
tic). Given the temperature of background gravitons, their en- 
tropy can be calculated as 

_2^_fc4 

■^grav " 



45 c^r?^ " 



(21) 



where ggi-av = 2. 

Figure|2]shows as a function of temperature. The func- 
tion is well known for temperatures below about lO'^eV, 
but is not known at higher temperatures. Previous esti- 
mates of the background graviton entropy have assumed 
g..s(fni.nA) g.smyV) = 106.75 (IFrampton et alJ 120091: 
IFrampton & Kephar^ 120081) . but this should be taken as a 
lower bound on §*s(fpianck) yielding an upper bound on Tgrav 

and ^grav ■ 

To get a better idea of the range of possible graviton tem- 
peratures and entropies, we have adopted three values for 
§*s(fpianck)- As a minimum likely value we use g^,s = 200 
(Figure |2l thick blue line), which includes the minimal set 
of additional particles suggested by supersymmetry. As our 
middle value we use = 350, corresponding to the linear 
extrapolation of g^,s in log(r) to the Planck scale (Figure |2l 
gray line). And as a maximum likely value we use = 10^, 
corresponding to an exponential extrapolation (Figure|2] thin 
blue line). 

The corresponding graviton temperatures today are (Equa- 
tion [SO]) 

7g,av = 0.6l!!]j^K. (22) 

Inserting this into Equation (|2T]) we find the entropy in the 
relic graviton background to be 



^+0.2 _o 

Vav = 1.7 X 10 km , 
5„„v = 6.2x 10^^-" k. 



(23) 
(24) 

It is interesting to note the possibility of applying Equa- 
tion ( |20] | in reverse, i.e., calculating the number of relativistic 
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10'" 10'° 10" 

temperature [eV] 

Fig. 2. — Number of relativistic degrees of freedom g^s as a function 
of temperature, computed using the prescription given bv lColeman & Rooi 
120031) . All the particles of the standard model are relativistic at T > lO'^ eV 

and g,5(I0'- eV) = 106.75. The value of is not known above T ~ lO'-. 
To estimate plausible ranges of values, we extrapolate g^s linearly (gray line) 
and exponentially (thin blue line) in log(r). The minimum contribution to 
gfS from supersymmetric partners is shown (blue bar) and taken to indicate 
a minimum likely value of g^s ^t higher temperatures (thick blue line). 



degrees of freedom at the Planck time using future measure- 
ments of the graviton background temperature. 

2.5. Dark Matter 

The most compelling interpretation of dark matter is as a 
weakly interacting superpartner (or weakly interacting mas- 
sive particle, WIMP). According to this idea, dark matter par- 
ticles decoupled from the radiation background at some en- 
ergy above the particle mass. 

If this interpretation is correct, the fraction of relativistic 
background entropy in dark matter at the time dark matter 
decoupled fdm dec is determined by the fraction of relativistic 
degrees of freedom that were associated with dark matter at 
that time (see Equation[T4li. 



Sdn 



S^S dm(^dm dec) 

-dm(^dm dec) 



■^non-dm rad 



(25) 



t *5 non- 



This can be evaluated at dark matter decoupling, or any time 
thereafter, since both s^,^ and inon-dm rad are adiabatic (oc a~^). 

We are unaware of any constraint on the number of su- 
perpartners that may collectively constitute dark matter. The 
requirements that they are only weakly interacting, and that 
they decouple at a temperature above their mass, are proba- 
bly only satisfied by a few (even one) species. Based on these 
arguments, we assume g^s dmikm dec) < 20 and g*s(fdm dec) > 



106.75 which yields the upper limit 

dm(^dm dec) 



dm dec 



< 



(26) 



On the other hand there may be many more degrees of free- 
dom than suggested by minimal supersymmetry. By extrap- 
olating g^,s exponentially beyond supersymmetric scales (to 
10^^ eV), we find §*s(fdmdec) < 800. In the simplest case, 
dark matter is a single scalar particle so g^s dm(fdm dec) > 1 
and we take as a lower limit 



g*S dm(fdm dec) 



> 



1 



6*o non- dm(fdm dec) ^ 800 

Inserting this into Equation dZST l at the present day gives 



(27) 



^dn 



= 5 X 10^^' /ti 



(28) 



where we have used the estimated limits given in Equations 
( |26] | and dZTl ) and taken inon-dm rad to be the combined entropy 
of neutrinos and radiation today (Equations [T2l and [Tsl). The 
corresponding estimate for the total dark matter entropy in the 
observable universe is 



5dm = 2x 10' 



88±1 



(29) 



As with our calculated neutrino entropy, our estimates here 
carry the caveat that we have not considered changes in the 
dark matter entropy associated with gravitational structure 
formation. 

2.6. Stellar Black Holes 

In the top panel of Figure |3] we show the stellar initial mass 
function (IMF) parameterized by 



drii. 



dlog{M) 



oc 



M 



(30) 



with a = -1.35 at M <0 .5Mc^ and a = -2.35![jj^ at M > 
0.5Mq (lElmegreenll2007b . We also show the present distri- 
bution of main-sequence stars, which is proportional to the 
initial distribution for M < IMq, but which is reduced by a 

factor of (M/Mq)"^-^ for heavier stars ( iFukugita & PeeblesI 
12001 . 



dn 



dn„. 



present 



dlogiM) 



dlog{M) 
rflog(M) 



(jl) 



-2.5 



for M < IM, 
forM> IM, 







(31) 







The initial and present distributions are normalized using 
the present cosmic densit y of stars, fi* = 0.0027 ± 0.0005 
(IFukugita & Peeblesll200l . 

The yellow fill in the top panel represents stars of mass 
IMq <M< SMq, which died leaving white dwarf remnants 
of mass M < IAMq (yellow fill, bottom panel). The 
blue fill represents stars of mass 8Mq <M <25Mq, which 
died and left neutron star remnants of mass 1 AMq <M< 
2.5Mq. The light gray area represents stars of mass 25M0 < 
M < 42Mq which became black holes of mass 2.5Mq <M < 
15M0 via super novae (here we use the sim plistic final-initial 
mass function of lFrver & Kalogeral ( 1200 lb ). Stars larger than 
~ 42Mq collapse directly to black holes, without supernovae, 
and therefore retain most of their mass ( dark gray regions; 
iFrver & Kalogera.200h.Heger et alJlIOOSi) . 
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progenitors 




remnants 



3 ixio" k- ^ 



10-- 

10^ - 
0.1 



BHs via 
NSs fallback 



BHs via 
direct 
collapse 

100 



- 10"' 



Fig. 3. — Progenitors in the IMF (top panel) evolve into the distribution of 
remnants in the bottom panel. The shape of the present main-sequence mass 
function differs from that of the IMF (top panel) by the stars that have died 
leaving white dwarfs (yellow), neutron stars (blue), and black holes (light 
and dark gray). The present distribution of remnants is shown in the bottom 
panel. Black holes in the range 2.5Mq <^M <^ 15Mq (light gray) have been 
observationally confirmed. They form from progenitors in the range 25Mq <^ 
M ~ 42Mq via core collapse supernova and fallback, and we calculate their 

q-j+0.6 

entropy to be 5.9 X 10 -i ^i. Progenitors above about 42 Mq may evolve 
directly to black holes without significant loss of mass (dark gray) and may 
carry much more entropy, but this population has not been observed. The 
green curve, whose axis is on the right, shows the mass distribution of stellar 
black hole entropies in the observable universe. 



5sMBH = 10"" - 10""/t. Bel ow we use the SMBH mass func- 
tion as measured recently bv lGraham et alj ( 120071) . Assuming 
a three-parameter Schechter function 



dn 



dlog{M) 



Q+l 



T7- exp 



M 

M, 



1-^ 



M 



(34) 



(number density per logarithmic mass interval) they find 0* = 
0.0016 ±0.0004 M/7c"^ = 2.9 ±0.7 x 10^ Mq, and a = 
-0.30 ± 0.04. The data and best-fit model are shown in black 
in Figure m 



10" 



10- - 




Fig. 4. — The black curve, whose axis is on the left, is the SMBH mass 
function from Graham et al. (20()'3), i.e., the number of supermassive black 
holes per Mpc^ per logarithmic mass interval. The green curve, whose axis is 
on the right, shows the mass distribution of SMBH entropies in the observable 
universe. 



Integrating Equation (|3]l over stellar black holes in the range 
M < 15Mq (the light gray fill in the bottom panel of Figure 
O we find 

1 7+0.(i _T 

ssBH (m<i5Mq)= 1-6 X 10 -' ^-km , (32) 
SsBH (m<i5Mq) = 5.9 X lO"'--" k, (33) 

which is comparable to previous estimates of the stellar black 
hole entropy (see Table [T]!. Our uncertainty is dominated by 
uncertainty in the slope of the IMF, but also includes uncer- 
tainty in the normalization of the mass functions and uncer- 
tainty in the volume of the observable universe. 

If the IMF extends beyond M > 42M0 as in Figure [3] then 
these higher mass black holes (the dark gray fill in the bot- 
tom panel of Figure O may contain more entropy than black 
holes of mass M < 15 Mq (Equation |32]i. For example, if the 
Salpeter IMF is reliable to M = 140 Mq (the Eddington limit 
and the edge of Figure |3]l, then black holes in the mass range 

42- 140 Mq would contribute about 3.1 x lO''*''''*'' k to the en- 
tropy of the observable universe. Significantly less is known 
about this potential population, and should be considered a 
tentative contribution in Table [T] 

2.7. Supermassive Black Holes 

Previous estimates of the SMBH entropy jPenrosd 120041 : 
iFrampton et"aLl 120091: iFrampton & Kepharli l2008h have as- 
sumed a typical SMBH mass and a number density and yield 



We calculate the SMBH entropy density by integrating 
Equation (O over the SMBH mass function, 

s= — — / M^ — d\og{M)- (35) 

cTi J \d\og{M)) ' 

The integrand is plotted using a green line in Figure |4] show- 
ing that the contributions to SMBH entropy are primarily due 
to black holes around ^ 10''Mq. The SMBH entropy is found 
to be 

SSMBH = 8.4!12^ 1023^m-\ (36) 
5sMBH = 3.l!?:?xlO'°4^. (37) 

The uncertainty here includes uncertainties in the SMBH 
mass function and uncertainties in the volume of the observ- 
able universe. This is at least an order of magnitude larger 
than previous estimates (see Table [T]i. The reason for the dif- 
ference is that the (Graham et al. 2007) SMBH mass func- 
tion contains larger black holes than assumed in previous es- 
ti mates. 

iFramptonI (l2009allbl) has suggested that intermediate mass 
black holes in galactic halos may contain more entropy than 
SMBHs in galactic cores. For example, according to the mas- 
sive astrophysical compact halo object (MACHO) explana- 
tion of dark matter, intermediate mass black holes in the mass 
range 10^ - 10^ Mq may constitute dark matter. Assuming 
10^ Mq black holes, these objects would contribute up to 
lO'"^ k to the entropy of the observable universe (IFramptonI 
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Fig. 5. — Whether or not the total black hole entropy is dominated by 
SMBHs depends on the yet-unquantified number of intermediate mass black 
holes. 

I2009bl) . Whether or not this is so depends on the number den- 
sity and mass distribution of this population. Figure |5] com- 
bines Figures|3]and|4]and shows what intermediate black hole 
number densities would be required. 

3. THE ENTROPY OF THE COSMIC EVENT HORIZON AND ITS 
INTERIOR 

In this section we calculate th e en tropy budget for scheme 
2 (refer to discussion in Section [LTT i. Scheme 2 differs from 
scheme 1 in two ways; first, along with the components pre- 
viously considered (and listed in Table [TJ, here we consider 
the CEH as an additional entropy component; and second, the 
volume of interest is that within the event horizon not the par- 
ticle horizon (or observable universe). 

The proper distance to the CEH is generally time- 
dependent, increasing when the universe is dominated by an 
energy component with an equation of state w > -I (radia- 
tion and matter) and remaining constant when the universe 
is dark energy dominated (assuming a cosmological constant, 
w = -1). Since our universe is presently entering dark energy 
domination, the growth of the event horizon has slowed, and 
it is almost as large now as it will ever become (bottom panel 
of Figure[T]i. In the Appendix, we calculate the present radius 
and volume of the CEH 

7?CEH = 15.7 ± 0.4 GZyr, (38) 

VcEH = 1 -62 ± 0. 12 X lO'* G;>'r^ 

= 1.37 ±0.10 X 10™ (39) 

We also calculate the prese nt entropy of the CEH (following 
iGibbons & Hawkingl[l977h . 

kc^A 
kc^ 

= ^^^CEH (40) 

= 2.6 ±0.3 X 10'^^ k. 

Entropies of the various components within the CEH are 
calculated using the entropy densities i, from Section|2] 

Si = SiVcEH (41) 



Table |2] shows that the cosmic event horizon contributes al- 
most 20 orders of magnitude more entropy than the next 
largest contributor, supermassive black holes. 

4. DISCUSSION 

The second law of thermodynamics holds that the entropy 
of an isolated system increases or remains constant, but does 
not decrease. This has been applied to the large-scale universe 
in at least two ways (Equation [T] and |2li. The first scheme 
requires the entropy in a comoving volume of the universe 
to not decrease. The second scheme requires the entropy of 
matter contained within the event horizon, plus the entropy of 
the event horizon, to not decrease. 

We have calculated improved estimates of the current en- 
tropy budget under scheme 1 (normalized to the current ob- 
servable universe) and scheme 2. These are given in Tables[T] 
andlH respectively. 

The entropy of dark matter has not been calculated previ- 
ously. We find that dark matter contributes 10*^**^' k to the 
entropy of the observable universe. We note that the neu- 
trino and dark matter estimates do not include an increase 
due to their infall into gravitational potentials during struc- 
ture formation. It is not clear to us a priori whether this non- 
inclusion is significant, but it may be since both components 
are presently non-relativistic. This should be investigated in 
future work. 

Previous estimates of the relic graviton entropy have as- 
sumed that only the known particles participate in the rel- 
ativistic fluid of the early universe at f > fpianck- In terms 
of the number of relativistic degrees of freedom, this means 
g^,s — > 106.75 at high temperatures. However, additional par- 
ticles are expected to exist, and thus g^,s is expected to become 
larger as f — s- fpianck- In the present work, we have calculated 
the relic graviton entropy corresponding to three high-energy 
extrapolations of (constant, linear growth and exponential 
growth) and reported the corresponding graviton temperatures 
and entropies. 

In this paper, we have computed the entropy budget of the 
observable universe today Sobsit = to)- Figure |6] illustrates the 
evolution of the entropy budget under scheme 1, i.e., the en- 
tropy in a comoving volume (normalized to the current ob- 
servable universe). For simplicity, we have included only the 
most important components. 

At the far-left of the figure, we show a brief period of in- 
flation. During this peri od all of the energy is in the inflaton 
( lGuthlI981l:lLindelI982h . whic h has very few degrees of free - 
dom and low entropy (blue fill: lLmdel2009l:[Steinhardtl2009l) . 
Inflation ends with a period of reheating somewhere between 
the Planck scale (IQ-'^^s) and the GUT scale (10"^^i), dur- 
ing which the inflaton's energy is transferred into a relativistic 
fluid (yellow fill). During reheating, the entropy increases by 
many orders of magnitude. After reheating, the constitution 
of the relativistic fluid continues to change, but the changes 
occur reversibly and do not increase the entropy. 

After a few hundred million years (^ lO^'i), the first stars 
form from collapsing clouds of neutral hydrogen and helium. 
Shortly thereafter the first black holes form. The entropy in 
stellar black holes (light gray) and SMBHs (dark gray) in- 
creases rapidly during galactic evolution. The budget given 
in Table [T] is a snapshot of the entropies at the present time 
(4.3 X lO'^i). Over the next IQ-^s, the growth of structures 
larger than about 10'"* Mq will be halted by the acceleration 
of the universe. Galaxies within superclusters will merge and 
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TABLE 2 

Entropy of the Event Horizon and the 
Matter Within it (Scheme 2 Entropy 
Budget) 



Component 



Entropy 5 [k] 



Cosmic Event Horizon 
SMBHs 

Stellar BHs (2.5 - 15 Mq) 

Photons 

Relic Neutrinos 

WIMP Dark Matter 

Relic Gravitons 

ISM and IGM 

Stars 



Total 



2.6 ±0.3 X 10122 
1.2+11 X 10103 

2.2 X I0<:2 

2.03 ±0.15 X 10**!* 
1.93 ±0.15 X 10**!* 
6x 10»«±i ^ 

2.3 X io'*'^^:i 

2.7 ±2.1 X 10*1 

3.5 ± 1.7 X 10''^ 

2.6 ±0.3 X 10'2^ 



Tentative Components: 
Massive Halo BHs (10^ Mq) 
Stellar BHs (42-14OM0) 



10104 
1.2 X 10' 



scalefactor 
1 10'°'° 10" 



Note. — This budget is dominated by the cos- 
mic event horizon entropy. While the CEH entropy 
should be considered as an additional component 
in sche me 2, it also c orresponds to the holographic 
bound I't Hooftll 19931) on the possible entropy of the 
other components and may represent a significant 
overestimate. Massive halo black holes at 10^ Mq 
and stellar black holes in the range 42—140 Mq are 
included tentatively since their existence is specula- 
tive. 




lO"" 10'" 10°' 
lime [s] 

Fig. 6. — The entropy in a comoving volume (normalized to the present ob- 
servable universe). This figure illustrates the time-dependence of the scheme 



scalefactor 
10-'° 10° 




c 10''° - 



1 entropy budget. N.B. 10" 



: 1 googolplex. 



lime [s] 



Fig. 7. — Entropy of matter within the CEH, and the entropy of the cosmic 
event horizon. This figure illustrates the time dependence of the scheme 2 
entiopy budget. Note: the horizontal axis is shorter than in Figure[6] 



those in the outer limits will be ejected. The final masses of 
SMBHs will be - IO'^Mq ( I Adams & LaughUnllT997l) with 
the entropy dominated by those with M IQ^^Mq. 

Stellar black holes will evaporate away into Hawking radi- 
ation in about 10^°s and SMBHs will follow in 10"°i. The 
decrease in black hole entropy is accompanied by a compen- 
sating increase in radiation entropy. The thick black line in 
Figure |6] represents the radiation entropy growing as black 
holes evaporate. The asymptotic future of the entropy bud- 
get, under scheme 1, will be radiation dominated. 

Figure |7] illustrates the evolution of the entropy budget un- 
der scheme 2, i.e., the entropy within the CEH, plus the en- 
tropy of the CEH. 

Whereas in scheme 1 we integrate over a constant comov- 



ing volume, here the relevant volume is the event horizon. 
The event horizon is discussed in some detail in the Ap- 
pendix. During radiation domination, the comoving radius of 
the CEH is approximately constant (the proper distance grows 
as ^CEH c>c fl) and in the dark energy dominated future, it is a 
constant proper distance {Rceh = constant). The few logarith- 
mic decades around the present time cannot be described well 
by either of these. 

Since the event horizon has been approximately comoving 
in the past, the left half of Figure |7] is almost the same as 
in Figure |6] except that we have included the event horizon 
entropy (green fill). The event horizon entropy dominates 
this budget from about lO''^^. 

After dark energy domination sets in, the CEH becomes 
a constant proper distance. The expansion of the universe 
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causes comoving objects to recede beyond the CEH. On aver- 
age, the number of galaxies, black holes, photons etc. within 
our CEH decreases as . The stellar and SMBH entropy 
contained within the CEH decreases accordingly (decreasing 
gray filled regions). 

The decreasing black hole entropy (as well as other compo- 
nents not shown) is compensated by the asymptotically grow- 
ing CEH entropy (demonst rated explicitly for a range of sce- 
narios in iDavis et aril2003h . a nd thus the second law of th er- 
modynamics is satisfied. See lEgan & LineweaveJ (120101 in 
preparation) for further discussion of the time-dependence of 
the entropy of the universe. 
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APPENDIX: THE OBSERVABLE UNIVERSE AND THE COSMIC 
EVENT HORIZON 

Here we calculate the radius and volume of the observable 
universe (for use in Section |2]i; and we calculate the radius, 
volume, and entropy of the CEH (for use in Section [3]l. We 
use numerical methods to track the propagation of errors from 
the cosmological parameters. 

The radius of the observable universe (or particle horizon) 



^obs = a(t) 



=0 a(t') 



-dt'. 



(42) 



Here a(t) is the time-dependent scalefactor of the universe 
given by the Friedmann equation for a flat cosmology 



da 
'dt 



fl- a 



(43) 



Hubble's constant and the matter density parameter are taken 
from lSeliak elall (l2006h : h = H/\00kms-^ Mpc"' =0.705 ± 
0.013, u}„, = n„h^ = 0.136± 0.003. The radiation density 
is calculated from the observed CMB temperature, Tcmb = 
2.725 ±0.002 K (iMaflier et al.lll999l) . using = i^f^- 
The vacuum energy density parameter is determined by flat- 
ness, r^A = 

A distribution of /?obs values is built up by repeatedly eval- 
uating Equation (l42l) at the present time (defined by fl(fo) =1) 
using cosmological parameters randomly selected from the 
allowed region of h - uj,„ - Tcmb parameter space (assuming 
uncorrected Gaussian errors in these parameters). We find 



7?obs = 46.9 ±0.4 Glyr 



(44) 



with an approximately Gaussian distribution. The quoted con- 
fidence interval here, and elsewhere in this Appendix, is la. 
The volume of the observable universe Vobs is calculated using 
the normal formula for the volume of a sphere. 



yobs = 43.2±1.2x lO'^G/yr 
= 3.65±0.10x 10' 



(45) 



See Figure [8] Uncertainty in /?obs and Vobs is predomi- 
nantly due to uncertainty in w,„ however h also makes a non- 
negligible contribution. 




5 1 



10 20 30 40 50 

volume of the particle horizon, V^^^ [10* Glyr^] 

Fig. 8. — Eight hundred realizations of Vqi,, ™d Vceh indicate the volume 
of the observable universe is 43.2 ± 1.2 X lO"* Glyr^ (horizontal axis) and the 
volume of the cosmic event horizon is Vceh = 1 -62 it 0. 12 X lO** Glyr^ (verti- 
cal axis). We note that there is only a weak correlation between uncertainties 
in the two volumes. 

The radius of the CEH at time t is given by integrating along 
a photon's world line from the time t to the infinite future. 



],=,„,„, a(t) 



(46) 



This integral is finite because the future of the universe is dark 
energy dominated. Using the same methods as for the observ- 
able universe, we find the present radius and volume of the 
CEH to be 



and 



7?CEH= 15.7 ±0.4 Glyr, 

VcEH=1.62±0.12x 10"^ Glyr\ 
= 1.37±0.10x 10™ ml 



(47) 



(48) 




0.68 0.70 
Hubble's constant, Hn [100 km 



0.7S 
Mpc" 



0.74 



Fig. 9. — We find 5ceh = 2.6 it 0.3 x lO'^^ k, in agreement with previous 
estimates 5ceh ~ lO'^^ k IBousso et alj|2007l) . Uncertainties in Sceh come 
from uncertainties in i?cEH. which are almost exclusively due to uncertainties 
in h. 



The entropy of the CEH is calculated using the 
Bekenstein-Hawking horizon entropy equation as suggested 
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Fig. 10. — Proper distance to the event horizon is shown as a function of 
time. The vertical gray line represents the present age of the universe (and its 
width, the uncertainty in the present age). During dark energy domination, 
the proper radius, proper volume, and entropy of the CEH will monotonically 
increase, asymptoting to a constant. 



Uncertainty in the CEH radius, volume, and entropy are dom- 
inated by uncertainties in Hubble's constant (Figure|9]l. 

The CEH monotonically increases, asymptoting to a con- 
stant radius and entropy slightly larger than its current value 
(see Figure [TOb. We calculate the asymptotic radius, volume, 
and entropy to be 

RcEHit (X)) = 16.4 ± 0.4 Glyr 

= 1.55 ± 0.04 X lO^Sii (50) 



VcEH(f oo) = 1 .84 ± 0. 15 X 10"* Glyr^ 

= 1.56±0.13 X 10™m^ (51) 



by lGibbons & Hawkingl(ll977h . 



kc^A kc^ , 



:2.6±0.3 X 10 



122 



(49) 



ScEHit -> oo) = 2. 



l±0.16x lO'-n. 



(52) 
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